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We report on some unusual behavior of the measured current- voltage characteristics 
(CVC) in artificially prepared two-dimensional unshunted array of overdamped Nb — 
AlO x — Nb Josephson junctions. The obtained nonlinear CVC are found to exhibit 
a pronounced (and practically temperature independent) crossover at some current 



I cr = ^2^: — lj Ic from a resistance R dominated state with Vr = Ry I 2 — Ic below 
I cr to a capacitance C dominated state with Vc = \J ^qVI ~ Ic above I cr . The 
origin of the observed behavior is discussed within a single-plaquette approximation 
assuming the conventional RSJ model with a finite capacitance and the Ambegaokar- 
Baratoff relation for the critical current of the single junction. 

PACS numbers: 74.50.+r, 74.81.Fa 

Among many different properties which ular mechanism and collateral effects (such as 
can be studied using highly ordered two- finite size of the single junction and/or the ar- 
dimensional arrays of Josephson junctions ray, thermal fluctuations, quasi-particle con- 
probably one of the most interesting (and tributions, etc 4r --). At the same time, recall 
important for their potential applications) is that only sufficiently overdamped Josephson 
their transport behavior, reflecting manifes- junctions (with nonhysteretic CVC) and their 
tation of numerous dissipation mechanisms in arrays can be effectively used in rapid sin- 
the arrays (see, e.g.,- - - and further references gle flux quantum (RSFQ) logic circuits and 
therein) via nonlinear current- voltage charac- programmable Josephson voltage standards 
teristics (CVC) of the general form V oc [I — (see, e.g.,— ~— and further references therein). 
IciT)] 11 ^ with a power exponent a ranging In this paper we report our results on CVC 
from a < 1 to a > 1 depending on the partic- for SIS type array of strongly overdamped 
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Nb — A10 X — Nb junctions at different tem- 
peratures. We observed quite a pronounced 
crossover at some current L. r = (77k 1)1, 



Lcr \2/3 c V c 

between a resistance R dominated state (be- 
low 7 cr ) and a capacitance C dominated state 
(above I cr ) which could be utilized as a versa- 
tile R-C switch within Josephson electronics. 
High quality ordered SIS type unshunted ar- 
ray of overdamped Nb — AlO x — Nb junctions 
has been prepared by using a standard pho- 
tolithography and sputtering technique 1 -. It 
is formed by loops of niobium islands linked 
through 100 x 150 tunnel junctions. The 
unit cell of the array has square geometry 
with lattice spacing a ~ 46fim and a sin- 
gle junction of 5 x 5/im 2 . The criti- 
cal current for the junctions forming the ar- 
rays is Ic{T) = 150 fi A at T = 1.7K. Given 
the values of the junction quasi-particle resis- 
tance R = 10Q and capacitance C = 1.2 fF, 
the circuit frequency and dissipation measur- 
ing Stewart-McCumber parameter are esti- 
mated to be uj rc = 1/CR ~ 10 14 Hz and 
P C (T) = 27rCJ f o Jc(r) ~ 0.05 at T = 1.7K, re- 
spectively. The parameters of the array are as 
follows, I C a(1-7K) = 1.2mA and R A = 25Q. 
The measurements were made using home- 
made experimental technique with a high- 
precision nanovoltmeter-. The temperature 
dependence of the normalized critical cur- 
rent of the array Ica(T)/Ica(0) is shown 
in Fig.l. Observe that it is very well fit- 
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FIG. 1. (Color online) The temperature depen- 
dence of the critical current for the array of over- 
damped Nb — AlO x — Nb Josephson junctions. 
The solid line is the best fit using the single- 
junction Ambegaokar-BaratofF relation. 

ted (solid line) to the Ambegaokar-Baratoff 
relation 14 for the critical current of a sin- 
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gle junction I c{T) = lc(0) 
where A(T) = A (0) tanh (2.2 ^f^) is the 
analytical approximation- of the BCS gap 
parameter (valid for all temperatures) with 
A(0) = 1.76k B T c ^ 1.5meV and J c (0) = 
7rA(0)/2ei? ~ 150/iA This remarkable ex- 
perimental fact suggests a rather strong coor- 
dinated response from all the junctions form- 
ing the array and allows us to substantially 
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simplify the analysis of the obtained results phase difference across the ith junction, and 

by considering the properties of a single junc- $o is the magnetic flux quantum. For any 

tion (or plaquette, see below). For this pur- particular solution <f>i(t) of this equation at 

pose, the initial CVC data for array were 1^0, the resulting CVC of the RSJ model 

rescaled by directly introducing the critical is given by the time average of the voltage 

current of the single junction Ic- Some typ- (r = 2-n juj is properly defined period) V(I) = 

ical results of the normalized rescaled CVC ^ J T ' dtU^f). For the resistance dominated 

taken at different temperatures are shown in situation (when the capacitance related ef- 

Fig.2. The solid lines through the data points fects can be totally neglected, that is when 

are the best fits according to the following urc t ^ 1) ; the RSJ model has a well-known 

expressions: V R = V Q ^J(P - for I c < solution 9 &(t) = 2 tan" 1 tan Of) - ^ 



/ < I cr and V c = V y(I — Ic)/h for / > I cr with u = I 2 — I c , which brings about 

with V = SOmV and J = 3mA It is inter- V R = ^ Jo dt (^) = Ry/P - I c for R- 

esting to point out that, according to Fig.2, dominated CVC. As we can see, this depen- 

the crossover shows almost a universal behav- dence exactly corresponds to the fitting ex- 

ior taking place around [I — Ic(T)]/Iq ~ 0.4 pression for the observed CVC below I cr as- 

for all temperatures. To understand the ob- suming the Ohmic relation Vq = RIq. Let 

served behavior of the CVC in our array, in us turn now to the opposite situation and 

principle one would need to analyze in detail consider the capacitance dominated regime 

the dynamics of the array. However, as we when the resistance related effects can be to- 

have previously reported^"—, because of the tally neglected (that is when ujrct 1). 

well-defined periodic structure of our array, it In this case, the first integral of the RSJ 



is reasonable to expect that our experimental model reads ^ = y/2u p <J + j^<t>i + cos 4>i 
results can be quite satisfactory explained by where u p = J2eIc/hC is the plasmon fre- 



analyzing the dynamics of a single unit cell quency, and C\ is the integration constant, 

of the array. In our calculations, the unit Unfortunately, for 1^0, this equation can 

cell is a plaquette containing four identical not be solved exactly. Hence, let us con- 

Josephson junctions. By analogy with the sider its approximate solution assuming that 

resistively shunted junction (RSJ) model 9 , 4>i{t) = | + 6{t) with 6{t) |. By fix- 

the total current in the plaquette reads^ ing the arbitrary constant as C\ = — |J, 

I = I C {T) sin0i(t) + + Q ir^-- Here within this approximation, we obtain 6(t) = 

4>i(t) is the gauge- invariant superconducting VtH 2 with Vt = tu p ^J (I — Ic)/2I C , which in 
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FIG. 2. (Color online) The normalized rescaled current-voltage characteristics for unshunted array 
of overdamped Nb — AlO x — Nb Josephson junctions taken at various temperatures along with the 
best fits (solid lines) using expressions for Vr and Vq (see text). Here, Iq = 3mA and Vq = 30mV. 
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turn results in the following explicit form of 
C-dominated CVC, V c = ^ Jo dt (f) = 
yjl — Iq- As we can see, this depen- 
dence exactly corresponds to the fitting ex- 
pression for the observed CVC above I cr as- 
suming V = u p <\j2eIo/hC = huRc/2e and 
I = Vq/R— Iq/ (3c for the normalization pa- 
rameters. Furthermore, given the experimen- 
tal values for R and C, we obtain Vq ~ 30mV 
and Iq ~ 3mA. In turn, from the obvious 
identity Vn(I cr ) = Vc(I cr ), we readily obtain 
I cr = — l)Ic for the crossover current. 
Finally, a comment is in order regarding the 
employed here simplified model of 2D array. 
It should be noted that we completely ig- 
nored inductance (geometry) related effects 
which seem to be of less importance for the 
interpretation of the observed crossover than 
dissipation induced factors (resistance and 
capacitance). However, for more adequate 
description of the flux dynamics in truly 2D 
systems, these effects should be taken into ac- 
count. Indeed, as accurate numerical simula- 
tions have revealed^i^, both self-inductance 
and mutual inductance effects will have a sig- 
nificant impact on the array's dynamic prop- 
erties through creation of rather strong self- 
induced magnetic fields. 

We are indebted to R.S. Newrock, P. Bar- 
bara and C.J. Lobb for helpful discussions 
and for providing SIS samples. This work 
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